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We use multiple scattering theory to evaluate the Q 2 dependence of the forward-backward 
asymmetry, A z , of the missing momentum distribution in quasielstic A(e, e'p) scattering, 
which is expected to be affected by color transparency. The novel features of our analysis 
are a consistent treatment of the structure of the ejectile state formed after absorption 
of the virtual photon by the struck proton and a careful evaluation of the background 
asymmetry induced by the nonvanishing real part of the proton-nucleon scattering am- 
plitude. We find that the absolute magnitude of A z is dominated by this background at 
the Q 2, s attainable at the Continuous Electron Beam Accelerator Facility. However, the 
Q 2 dependence of the asymmetry is sensitive to the onset of color transparency, whose 
observation in high statistics experiments appears to be feasible, particularly using light 
targets such as 12 C or 16 0. 
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The recent A(e, e'p) experiment carried out at SLAC by the NE18 collaboration M 
did not observe any substantial Q 2 dependence of the integrated nuclear transparency 
up to Q 2 < 7 GeV 2 . On the other hand, theory predicts that final state interactions 
(FSI) should vanish at asymptotically high Q 2 on account of color transparency (CT) 
0. |3|. From the point of view of multiple scattering theory (MST), CT corresponds to a 
cancellation between the rescattering amplitudes with elastic (diagonal) and inelastic (off- 
diagonal) intermediate states. The NE18 finding was anticipated in ref. @], where it was 
argued that in the integrated nuclear transparency the effect of the inelastic rescatterings 
of the struck proton is still weak in the region Q 2 < 7 GeV 2 , covered by NE18 and 
achievable by the forthcoming experiments at CEBAF. The cancellation of the diagonal 
and off-diagonal rescatterings depends upon the missing momentum [H || 0], which can 
be used to enhance the contribution from inelastic rescatterings and make the onset 
of CT observable at CEBAF. The enhancement and/or suppression of the contribution 
of inelastic rescatterings result in a non-vanishing forward-backward asymmetry of the 
missing momentum distribution , w(p), defined as 

A z ( X ,y) = J ^ TW _, (1) 

where N + (AL) is the number of events in the kinematical region x < p z < y {—y < p z < 
—x) and the z-axis is parallel to the (e, e') momentum transfer q. 

Unfortunately, besides CT, there is another competing source of asymmetry associated 
with FSI: the nonvanishing ratio between the real and the imaginary part of the elastic 
proton-nucleon (pN) scattering amplitude, a P N- The forward-backward asymmetry, as- 
sociated with a P N, disregarded in refs.|| || [?]], was recently found to be sizeable 0. The 
role played by the real part of the rescattering amplitude in electron nucleus scattering 
at high Q 2 has been also pointed out in ref. ||, within the context of an analysis of the 
inclusive data from SLAC. 

In this paper we give a quantitative evaluation of how much the CT signal can be 
obscured by the effect of a nonzero a P N- From the technical point of view, one must 
calculate A z taking into account the coupled- channel treatment of the off-diagonal reseat- 



terings and the nonzero a P N on the same footing. Besides including the effect of a P N, 
we improve upon the previous works treating the coherency properties of off-diagonal 
rescatterings in a more accurate manner. 

The missing momentum distribution equals 



f r i i 

w(p) = j^y J d i r 1 d i r l d , 'r2 ■ ■ ■ d 3 f A p P (fi, r 2 , • ■ ■ , r A \ r[, r 2 , ■ ■ • , r A ) 

xS*(f 1 ,f 2 , ■ ■ ■ ,f A )S(f 1 ,r 2 , ■ ■ ■ ,r A )exp[ip(f[ - fi)] , (2) 



where p p {r\, r 2 , • • • , f A ; f x , r 2 , • • • , f A ) is the proton A-body semidiagonal density matrix 
| T0| and 

S{r 1 ,r 2 ,---,r A ) = tttt x • (3) 

(p\E) 

In eq.@, ^(rj, f 2 , r A ) is the operator describing the evolution of the ejectile state 
during its propagation through the nuclear medium and the ejectile state \E) is a three- 
quark state formed from the struck proton after absorption of the virtual photon at 
position fi. In terms of the electromagnetic current operator J em , \E) = J em {Q)\p) = 
J2 \i) (i\Jem(Q)\p) — G ip (Q)\i) , where the sum includes the complete set of three-quark 

i i 

states and Gi P (Q) = (i\J e m(Q)\p) = (i\E) denotes both the electromagnetic form factor 
of the proton (|z) = \p)) and all the transition form factors associated with the electroex- 
citation processes e + p — > e' + i |. 

We will consider the forward-backward asymmetry at moderate missing momentum 
p < 200 MeV, where the main effect of short range NN correlations in the initial state 
is an overall renormalization of the momentum distribution |TT|, which does not affect 



the calculation of A z . Although iViV correlations also contribute to the integrated trans- 
parency |12], 113] , their effect is expected to be less important in A z . Neglecting exchange 
and dynamical correlations between the struck nucleon and the spectator nucleons, the 
proton z/-body semidiagonal density matrix, p(fi, r 2 , ■ ■ ■ , F„; fi, r 2 , ■ ■ ■ , r„) takes the fac- 
torized form p p (fi, r 1 )p(f 2 , ■ ■ ■ , r„), where p{r 2 , ■ • • , r u ) is the (y — l)-body diagonal density 
matrix which can be recast in the form |14|] 11^=2 /°(^)p(^2) ' ' ' i ^v) where p(f*i) = n A [fi)j A 



is the one-body density and the [y— l)-body distribution function p0| , g(f 2 , ■ • ■ , f v ), brings 



in the effect of NN correlations. We also neglect NN correlations among the spectator 
nucleons in this first study of longitudinal asymmetry. This implies that g(f 2 , • • • , r u ) = 1 
and one ends up with the results of ref. [8] that at moderate missing momenta p < 200 
MeV/c, the integrand of @ - the FSI-modified one-body density matrix - can be approx- 
imated as 

, , A 

P P (n,K) / d3ff 2 • ■ -d 3 r A Y[p{ri)S*(f[,f 2 , ■ ■ ■ , f A )S(f 1 ,f 2 , • ■ ■ , f A ) = p p (fi, r-[)$(ri, f x ) 
J j=i 

(4) 

where, Pp{ri,r'\) = \ 2 n 4 > n( r i) ( f > n('>i) is the proton shell model one-body density matrix 
and 4> n are the shell model wavefunctions. 

In the MST, the FSI operator S^S^g can be very schematically represented as S^S^ = 
1X^=2(1 — r*^,-)(l — fy), where the profile function f is an operator acting on the state 
describing the internal structure of the three-quark system [8] . Then, the calculation of 
$(r*i, f[) involves a sum of diagrams, as shown in Fig. 1. Every dotted line attached to the 
straight-line trajectory originating from the point r\ (f[) denotes a profile function f\j 
(f 1 */ .), where the subscripts lj(l'j) stands for the position of the j-th spectator and its 
projection on the trajectory. The interaction between the two trajectories, due to terms 
TijT^j like the one shown in Fig. lb, makes the FSI factor $(r!,r{) a non factorizable 
function of f% and f 1 . It can be shown that the non factorizable part of $(r*i, f\) is only 
connected with the incoherent elastic rescatterings, which come into play only at p > 200 
MeV ||. For this reason, at p < 200 MeV, we can neglect FijT*, . terms in the FSI 
operator S 3q and obtain the factorizable <3>(ri,r{) = S* oh (fi)S co h{r[) and the following 
expression for w(p), 

(5) 
where 

Scohiri) = J \[p{r 5 )d rj . (6) 

Finally, using the Glauber form of Ss q , the coupled channel MST series for S co h(r) in terms 

00 

of the i/-fold off-diagonal rescatterings can be cast in the form S co h(r) = I] S^ oh (f) , where 
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d 3 ?(j) n (r) exp[-ipr\S coh (f) 



the first few terms read 



■ ( ' I E \ 2 
S C oh(b, z) = 1 - exp[ik hp z] V 1 / dz x n A {b, zi) exp[ifc pil 2i - -£(&, zi, 

+ exp[zfcj lP z] ^ ^ 2 77^- [ dzin A (b, z^ exp[ik i2h zi - ~t(b, zx, z)a ilh ] 
ili2 1 1 \P\&) { 1 

oo 

x y dz 2 n A (b, z 2 ) exp[ik pi2 (z 2 - z x ) - -t(b, z 2 , zi)a i2i2 } 



+ (7) 

where r = (b, z), a' ik = <Jik—Sik<Ju, the matrix a = 2 J d 2 bT(b) is connected with the forward 
diffraction scattering matrix f = ia (notice that Re an = cr tof (iiV)) and t(b,z 2 ,zi) = 
dzn A (b, z) is the partial optical thickness. Eq.(|7]) shows contributions from up to two 
off-diagonal rescatterings and any number of elastic rescatterings in between, higher order 
off-diagonal rescatterings are implied by the dots. The onset of CT is controlled by the 
oscillating factors exp[iki n i n _ 1 z n ], taking into account the longitudinal momentum transfer 



associated with each off-diagonal transition i n -\N — > i n N [15] according to 
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rrij — mf . 

I l n l n — l 



where rrij denotes the mass of the proton excitation |z) and e is the energy of the struck 
proton in the laboratory frame. 

From the point of view of CT, the most important ingredient in the coupled-channel 
formalism is the Pomeron contribution to the matrix a. As in ref. [|J, we will construct 
a using the oscillator quark-diquark model of the proton. As a result, the Pomeron 
contribution to the matrix element can be written as 

Im f(kN -> iN) = Re af k = J dzd 2 p%{p, z)cr(p)* k (p, z) , (9) 

where are the oscillator wave functions describing the quark-diquark states and er(p) 
is the dipole cross section describing the interaction of the quark-diquark system with a 
nucleon, which was taken in the form a(p) = oq 1 — exp ( — "f?) • Following ref.||, we set 
do = 2a to t{pN) and adjust R Q to reproduce atotipN). For the diagonal iN — ► iN scatter- 
ing we take Re f(iN — > iN) — — Im an — | (a P p<7tot(pp) + Q-pn&totijpn)) ■ Even though this 



choice can be justified within the framework of the dual parton model we are fully 
aware that it should only be regarded as an estimate. The real parts Re f(iN — > kN) of 
the off-dagonal amplitudes are not known experimentally. We shall give an estimate of 
the sensitivity of our results to these quantities and the associated uncertainty. 

For the oscillator frequency of the quark-diquark system we use the value ui q r, = 0.35 
GeV, leading to a realistic mass spectrum of the proton excitations. Since we start with a 
a(p) which vanishes as p — > 0, our diffraction matrix satisfies the CT sum rules || |l7j] by 
construction, and leads to vanishing FSI in (e, e'p) at asymptotically large Q 2 . Moreover, 
with the above form of cr(p), we obtain a diffraction matrix yielding a realistic description 
of the mass spectrum observed in diffractive pN scattering [|17]]. For instance, the present 
model gives a value of the ratio between the diffractive and the elastic pN cross sections 
in agreement with the experimental data: &dif f(pp) / <?ei(pp) ~ 0.25 |18| . The effect of 
the truncation of the sum over the intermediate states has been evaluated by direct z- 
integration in eq. ([?]). This procedure is more accurate than the effective diffraction 
matrix approximation used in earlier works @, [|. The main difference is that the direct 
evaluation of the z-integrations in eq.([?J produces a somewhat weaker suppression of the 
contribution of heavy intermediate states and a somewhat faster onset of CT effects. 
Since our diffraction matrix gives a reasonable description of diffractive pN scattering, we 
believe that the multiple scattering approach employed in this paper is a good tool for a 
quantitative study of the onset of CT in (e, e'p) reactions. 

Besides the matrix a, the evaluation of w(p) requires the initial wave function \E). 
Assuming the dominance of the hard mechanism in the electromagnetic form factors [H| . 
it can be shown that the projection of the ejectile onto the subspace of hadronic mass 
eigenstates which satisfy the coherency constraint m* 2 — m 2 <Q 2 / 'iJ^m^, coincides with 
the similar projection of a compact state, with transverse size p ~ 1/Q [|20| , 17] . For 
this reason, the latter can be taken for \E), and in our analysis we use the initial wave 
function of the form (p\E) oc exp(— Cp 2 Q 2 ), with C = 1. It is worth noting that the 
missing momentum distribution is not sensitive to the specific choice of C as long as 
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C > 1/ Q 2 p 2 , where p Q denotes the position of the first node in the wave functions of the 
excited states satisfying the coherency requirement. In the region of 2 < Q 2 < 20 GeV 2 , 
that we discuss in the present paper, the numerical results are practically independent of 
the parameter C for C > 0.05. 

Before discussing our numerical results, a comment on the treatment of A z of refs. |5|, 
|7|, [Tl[] is in order. First, these authors have not accounted for the large contribution 
to A z coming from the nonzero a P N- Second, we disagree with the treatment of the 
initial state \E) suggested in ref.||. The formalism of the coupled channel eikonal wave 
equations used in H 0, ^TJ is equivalent to our coupled channel multiple scattering series of 
eq.(|7|). However, it has to be noticed that, within the Plane Wave Impulse Approximation 
(PWIA), eq. (0) combined with our definition of the ejectile wave function \E) gives 
da(A(e,e'i))/da(A(e,e'p)) = B 2 (p)/B 2 (p), where 

d 3 kn(k x , k y , k z )5 [k z - p z ^ G%{Q 2 ) , (10) 

and n(p) is the nucleon momentum distribution. Eq. ([10|) shows the standard kinematical 
correlation between the mass of the excited baryon state \i) and the missing momentum, 
which can produce either enhancement or suppression of the resonance abundance in the 
final state. The procedure used in refs.|5|, [7], ^1[ amounts to a redefinition of the ejectile 

state 

\E(p)) = E G *(G)!*j§|i> > (") 
leading to a p dependence in ejectile wave function ([TT|) that is not born out by the 
quantum-mechanical treatment of the Fermi motion. It corresponds to a double counting 
and evidently does not give the correct PWIA limit of the production cross sections. Con- 
sequently, the validity of the numerical results of refs. |^|, 0, ^T|] appears to be questionable 
even within the crude models employed to describe the diffraction matrix. 

The main results of our work are presented in Fig. 2. Following we have calculated 
A z (x, y) assuming parallel kinematics, i.e. taking in eq.(|l]) iV + = f% dp z w(p± = 0,p z ) and 
iV_ = fZy dp z w(p± = 0,p z ). Fig. 2 shows A z calculated for two kinematical windows, 



0, y) = (0, 200) and (x, y) = (50, 200) MeV, and two different targets, 16 and 40 Ca. The 
parameters of the pp- and pn-amplitudes were taken from the recent review of ref . |22| . 
The number of the included excited states and the off-diagonal rescatterings were equal to 
4 and 3 respectively. Contributions from higher excitations and rescatterings with v > 3 
have been found to be negligible at Q 2 < 20 GeV 2 . In the region of Q 2 < 5 GeV 2 the 
effect of CT is almost exhausted by the first excitation. To illustrate the contribution 
of the elasic intermediate state to A z , the results obtained using the Glauber model are 
also shown in Fig. 2. In this case the only source of the z-asymmetry of the missing 
momentum distribution is the nonzero a P N- As one can see from Fig. 2, in the region 
of Q 2 ~ 2 — 5 GeV 2 , relevant to the experimental study of CT effects at CEBAF, the 
inelastic intermediate states are only responsible for 25-30 % (for 16 0) and 10-15 % (for 
40 C a) of A z for both the kinematical windows in p z . The large contribution of the elastic 
rescatterings to A z can make it difficult to extract definite conclusions on the onset of 
CT from measurements of the magnitude of A z at CEBAF. According to our results, the 
situation is particularly unfavourable in the case of heavy targets. However, in light nuclei 
the Q 2 -dependence of A z is dominated by CT. 

Since CT effects on A z are small, assessing the posssibility of their observation at 
CEBAF energies requires a study of the sensitivity of the theoretical predictions to the 
unknown ratios between the real and the imaginary parts of the resonance-nucleon am- 
plitudes. At Q 2 < 10 GeV 2 , our calculations show a weak sensitivity of A z to variations 
of ctjjv for the diagonal iN — > iN rescatterings of the excited states i. However, the effect 
of the uncertainty associated with the a{iN —>■ fciV)'s is not negligible. The hatched area 
In Fig. 2 shows the band of variation of A z corresponding to variations of a(iN — > kN) 
between -0.5 and 0.5. It appears that, at least for nuclear mass number A ~ 10, the 
uncertainty of the theoretical predictions does not rule out the possibility of using the Q 2 
dependence of A z in the Q 2 range attainable at CEBAF to explore the onset of CT. 

In conclusion, we have performed a MST calculation of the longitudinal asymmetry A z 
in (e, e'p) scattering. Our analysis improves upon the previous works on this problem in 
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several aspects. We have for the first time taken into account the background asymmetry 
induced by the nonvanishing ratio between the real and the imaginary part of the proton- 
and resonance-nucleon amplitudes, which turns out to be very important and dominates 
the asymmetry at moderate Q 2 . Moreover, we have developed a consistent treatment 
of the ejectile state, which is also free of problems of double counting. The onset of 
CT, connected with the occurrence of inelastic rescatterings, has been studied using a 
realistic model of the diffraction scattering matrix, without making approximations in 
the treatment of coherency effects in the coupled channel multiple scattering series. Our 
results suggest that the Q 2 dependence of A z , if measured in high statistics experiments 
off light targets, can be used as a signature of the onset of CT in the kinematical range 
accessible at CEBAF. 
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Foundation and the INTAS grant 93-239. AAU acknowledges Prof. A. Zichichi and 
ICSC- World Laboratory for financial support. BGZ wishes to gratefully acknowledge the 
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Figure captions: 

Fig. 1 - The typical diagrams contributing to the FSI factor $(ri,r{) in the MST: (a) the 
diagram without the interaction between the two trajectories outgoing from fi and 
r[i (b) the diagram containing the interaction between the trajectories generated by 
the term f (6j)f *(6^-). 

Fig. 2 - The Q 2 -dependence of the longitudinal asymmetry A z (x,y) as given by Eq. (1) 
in the case of parallel kinematics p± = for 16 0(e, e'p) and 40 Ca(e,e'p) scattering. 
The solid lines show the results obtained with inclusion of the elastic and inelastic 
intermediate states within the coupled-channel MST, while the dashed lines were 
obtained in the one-channel Glauber model. The shaded areas show the limits of 
variations of Af ST for variations of a(iN — > kN) between -0.5 and 0.5. 
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